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Abstract. Neutrinoless double beta decay is a very impor- 
tant process both from the particle and nuclear physics point of 
view. From the elementary particle point of view it pops up in 
almost every model. In addition to the traditional mechanisms, 
like the light neutrino mass, A and rj terms etc we can have di- 
rect R-parity violating supersymmetric (SUSY) contributions. 
In any case its observation will severely constrain the existing 
models and will signal that the neutrinos are massive Majorana 
particles. From the nuclear physics point of view it is challeng- 
ing, because: 1) The relevant nuclei have complicated nuclear 
structure. 2) The energetically allowed transitions are exhaust 
a small part of all the strength. 3) One must cope with the 
short distance behavior of the transition operators, especially 
when the intermediate particles are heavy (eg in SUSY mod- 
els). Thus novel effects, like the double beta decay of pions in 
flight between nucleons, have to be considered. 4) The interme- 
diate momenta involved are about 100 MeV/c. Thus one has to 
take into account possible momentum dependent terms in the 
nucleon current. We find that, for the mass mechanism, such 
modifications of the nucleon current for light neutrinos reduce 
the nuclear matrix elements by about 25%, almost regardless 
of the nuclear model. In the case of heavy neutrino the effect 
is much larger and model dependent. Taking the above effects 
into account, the availabe nuclear matrix elements for the exper- 
imentally interesting nuclei A = 76, 82, 96, 100, 116, 128, 130, 
136 and 150 and the presently available experimental limits on 
the half-life of the 0z//3/3-decay we have extracted the following 
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new limits: {m^) < 0.3eV/c^ A'^ < 4.0 x 10~^ for the R-parity 
violating parameter with reasonable choices of the parameters 
of SUSY models, 



1. Introduction 

The nuclear double beta decay can occur whenever the ordinary (single) 
beta decay is forbidden due to energy conservation or greatly suppressed 
due to angular momentum mismatch. The exotic neutrinoless double beta 
decay {Qv (3 (3 — decay) is the most interesting since it violates lepton number 
by two units. It is a very old process. It was first considered by Furry 
[Q exactly half a century ago as soon it was realized that the neutrino 
might be a Majorana particle. It was continued with the work of Primakoff 
and Rosen |^ especially when it was recognized that kinematically it is 
favored by 10^ compared to its non exotic sister 2vj3(3-(iecay. When the 
corresponding level of the 10^^^ lifetime was reached and the process was 
not seen, it was tempting to interpret this as an indication that the neutrino 
was a Dirac particle. The interest in it was resurrected with the advent of 
gauge theories which favor Majorana neutrinos and through the pioneering 
work of Kotani and his group it was brought again to the attention 
of the nuclear physics community. To-day, fifty years later, Oi^/3/3-decay 
continues to be one of the most interesting processes. 

From a theoretical point of view it is the most likely, if not the only, 
process capable of deciding whether or not the neutrino is a Majorana 
particle, i.e. it coincides with its own antiparticle ||, ^, ^, ||, ||]. It is 
expected to occur whenever one has lepton number violating interactions. 
Lepton number, being a global quantity is not sacred, but it is expected to 
be broken at some level. In short this process pops up almost everywhere, 
in every theory. 

From a nuclear physics point of view calculating the relevant nuclear 
matrix elements it is indeed a challenge. First almost all nuclei, which can 
undergo double beta decay, are far from closed shells and some of them 
are even deformed. One thus faces a formidable task. Second the nuclear 
matrix elements represent a small fraction of a canonical value ( i.e. the 
matrix element to the energy non allowed transition double Gamow- Teller 
resonance or some appropriate sum rule). Thus effects which are normally 
negligible become important here. Third in many models the dominant 
mechanism for 0;//?/3-decay does not involve intermediate light neutrinos, 
but very heavy particles and one must be able to cope with the short 
distance behavior of the relevant operators and wave functions. 
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From the experimental point of view it also very challenging to measure 
the slowest perhaps process accessible to observation. Especially since it 
is realized that, even if one obtains only lower bounds on the life time for 
this decay, the extracted limits on the theoretical model parameters may 
be comparable, if not better, and complementary to those extracted from 
the most ambitious accelerator experiments. 

The recent superkamiokande results have given the first evidence of 
physics beyond the Standard Model (SM) and in particular they indicate 
that the neutrinos are massive particles. It is important to proceed further 
and find out whether the neutrinos are Dirac or Majorana particles. As 
we have mentioned there might be processes other than the conventional 
intermediate neutrino mechanism, which may dominate Oz//3/3-decay. It has, 
however, been known that whatever the lepton violating process is, which 
gives rise to this decay, it can be used to generate a Majorana mass for the 
neutrino jl^. The study of the Oi^/3/3-decay is further stimulated by the 
development of grand unified theories (GUT's) and supersymmetric models 
(SUSY) representing extensions of the SU{2)l(E>U{1) SM. The GUT's and 
SUSY offer a variety of mechanisms which allow the Oj^/3/?-decay to occur 

The best known possibility is via the exchange of a Majorana neutrino 
between the two decaying neutrons ||, ^, ^ |^, |l^, ^ |^. Nuclear physics 
dictates that we study the light and heavy neutrino components separately. 
In the presence of only left-handed currents for light intermediate neutrinos 
the obtained amplitude is proportional to a suitable average neutrino mass, 
which vanishes in the limit in which the neutrinos become Dirac particles. 
In the case of heavy Majorana neutrino components the amplitude is pro- 
portional to the average inverse neutrino mass, i.e. it is again suppressed. 
In the presence of right handed currents one has one can have a contribu- 
tion similar to the one above for heavy neutrinos but involving a different 
(larger) average inverse mass and some suppression due to the heaviness of 

Wr. 

It is also possible to have, in addition, interference between the leptonic 
left and right currents, jl — ju interference. In this case the amplitude in 
momentum space becomes proportional to the 4- momentum of the neutrino 
and, as a result, only the light neutrino components become important. 
One now has two possibilities. First the two hadronic currents have a 
chirality structure of the same kind Jl — Jr.- Then one can extract from 
the data a dimensionless parameter A, which is proportional to the square 
of the ratio of the masses of the L and R gauge bosons, n = {mL/mRf . 
Second the two hadronic currents are left-handed, which can happen via the 
mixing of the two bosons. The relevant lepton violating parameter 77 is now 
proportional to this mixing. Both of these parameters, however, involve the 
neutrino mixing and they are proportional to the mixing between the light 
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and heavy neutrinos. 

In gauge theories one has, of course, many more possibihties. Exotic 
intermediate scalars may mediate Oi'/J/J-decay These are not favored 
in current gauge theories and are not going to be further discussed. In 
superstring inspired models one may have siglet fermions in addition to 
the usual right handed neutrinos. Not much progress has been made on 
the phenomenological side of these models and they are not going to be 
discussed further. 

In recent years supersymmetric models are taken seriously and semire- 
alistic calculations are taking place. In standard calculations one invokes 
universality at the GUT scale, employing in all 5 parameters, and use 
the renormalization group equation to obtain all parameters (couplings 
and particle masses) at low energies. Hence, since such parameters are 
in principle calculable one can use 0:^/3/?-decay to constrain some of the 
R-parity violating couplings, which cannot be specified by the theory 
@ 1^, Recent review articles §, § give a detailed 

account of the latest developments in this field. 

From the above discussion it clear that one has to consider the case of 
heavy intermediate particles. One thus has to consider very short ranged 
operators in the presence of the nuclear repulsive core. If the interacting 
nucleons are point-like one gets negligible contributions. We know, however 
that the nucleons are not point like and they have structure described by a 
form factor, which can be calculated in the quark model or parameterized 
in a dipole shape. This approach was first considered by Vergados |Q 
adopted later by almost everybody. The resulting effective operator has a 
range somewhat less than the proton mass (see sect. 4 below). 

The other approach is to consider particles other than the nucleons 
present in the nuclear soup. For 0^ —^ 0+ the most important such particles 
are the pions. One thus may consider the double beta decay of pions in 
flight between nucleons, like 

TT^ — > TT^ e~ e~ , n > p Tr"*" e~ e 

This contribution was first considered by Vergados |Q and was found to 
yield results of the same order as the nucleon mode with the above recipe 
for treating the short range behavior. It was revived by the Tuebingen 
group l]!?! , [l8| in the context of R-parity violating interactions, in which it 
appears to dominate. 

The other recent development is the better description of nucleon cur- 
rent by including momentum dependent terms, such as the modification of 
the axial current due to PCAC and the inclusion of the weak magnetism 
terms. These contributions have been considered previously but 
only in connection with the extraction of the rj parameter mentioned above. 
Indeed these terms were very important in this case since they compete with 
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the p-wave lepton wave function, which, with the usual currents, provides 
the lowest non vanishing contribution. In the mass term, however, only 
s-lepton wave functions are relevant. So these terms have hitherto been 
neglected. 

It was recently found 1 23 that for light neutrinos the inclusion of these 
momentum dependent terms reduces the nuclear matrix element by about 
25%, independently of the nuclear model employed. For the heavy neutrino, 
however, the effect can be larger and depends on the nuclear wave functions. 
The reason for expecting them to be relevant is that the average momentum 
< g > of the exchanged neutrino is expected to be large . In the case of 
a light intermediate neutrino the mean nucleon-nucleon separation is about 
2 fm which implies that the average momentum < g > is about 100 MeV. 
In the case of a heavy neutrino exchange the mean intcrnuclcon distance is 
considerably smaller and the average momentum < g > is supposed to be 
considerably larger. 

Since Ov — j3[3 decay is a two step process, one should in principle con- 
struct and sum over all the intermediate nuclear steps, a formidable job 
indeed in the case of the Shell Model Calculations (SMC). Since, however, 
the average neutrino momentum is much larger compared to the nuclear ex- 
citations, one can invoke closure using some average excitation energy (this 
does not apply in the case of 2vj3(3 decays). Thus one need construct only 
the initial and final 0+ nuclear states. In Quasiparticle Random Phase Ap- 
proximation (QRPA) one must construct the intermediate states anyway. 
In any case it was explicitly shown, taking advantage of the momentum 
space formalism developed by Vergados that this approximation is 
very good ||2^, |2^. The same conclusion was reached independently by 
others [|8). 

Granted that one takes into account all the above ingredients in order 
to obtain quantitative answers for the lepton number violating parameters 
from the the results of Oz^/3/3-decay experiments, it is necessary to evalu- 
ate the relevant nuclear matrix elements with high reliability. The most 
extensively use methods are the SMC ( for a recent review see Q) and 
QRPA( for a recent review see |8j). The SMC is forced to use few single 
particle orbitals, while this restriction does not apply in the case of QRPA. 
The latter suffers , of course, from the approximations inherent in the RPA 
method. So a direct comparison between them is not possible. 

The SMC has a long history |2|, 1^, ||, |3|||, H in in double beta 
decay calculations. In recent years it has lead to large matrices calculations 
in traditional as well as Monte Carlo types of calculations |3^, ^ 

^ (For a more complete set of references see Ref. |^) and suitable 
effective interactions. 

There have been a number of QRPA calculations covering almost all 
nuclear targets H, H, ||, ||, ||, ||, ||, |o[ |l[ ||. We also have 
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seen some refinements of QRPA, like proton neutron pairing and inclusion 
of renormalization effects due to Pauli principle corrections [|5^, Q . 

The above schemes, in conjunction with the other improvements men- 
tioned above offer, some optimism in our efforts for obtaining nuclear ma- 
trix elements accurate enough to allow us to extract reliable values of the 
lepton violating parameters from the data. We will review this in the case 
of most of the nuclear targets of experimental interest (^^Ge, ^"^Se, ^^Zr, 
^°°Mo, ^^^Cd, i28Te, i30Te, ^^exe, iso^Vd). 



2. Theory 

2.1. Majorana neutrino mass mechanism 

We shall consider the QvPf3-Aecay process assuming that the effective beta 
decay Hamiltonian acquires the form: 

= ^ [e-7,(l - l,y,i] Jf + [e-7,(l + l.V^n] J'J + h.c. (2) 

where e and i/^j^ v^j^ are field operators representing the electron and the 
left handed and the right handed electron neutrinos in the weak interac- 
tion basis, respectively. We suppose that neutrino mixing does take place 
according to 

k=l k=l 

U^K,,+± u!:fN,„ (4) 
fc=i fc=i 

where, (Nk) are fields of light (heavy) Majorana neutrinos with masses 
ink {mk << 1 MeV) and Mk {Mk >> 1 GeV), respectively. The matrices 
and J/g^^'' are approximately unitary, while the matrices u'^'^ and 
LQ^ are very small (of order of the up quark divided by the heavy neu- 
trino mass scales) so that the overall matrix is unitary. ly'kjNk satisfy the 
Majorana condition: i^k^k = C V^, Nk^k — C N^. , where C denotes the 
charge conjugation and ^, S are phase factors (the eigenmasses are assumed 
positive) . 

We assume both outgoing electrons to be in the Si/2 state and consider 
only 0^ — > Oj transitions are allowed. For the ground state transition 
restricting ourselves to the mass mechanism we obtain for the Oi^/3/3-decay 
inverse half-life |, | |, 0, |, i, 

{T?/2]-' = Goi[|^^a4^^J> + ry^M.^^^^P + h^,^^'^'"'!^] (5) 
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The lepton-number non-conserving parameters, i.e. the effective neutrino 
mass < niij > and 77^ ,77^ are given as follows: 

< > =t iU'^^'f ^'^ - t (^eTV 2, ^, (6) 

1 1 

= (-' + ^')E (f^i^)' 2,^, (7) 

with TUp (me) being the proton (electron) mass, k is the mass squared ratio 
of Wl and Wr and e their mixing. Gqi is the integrated kinematical factor 
[H nuclear matrix elements associated with the exchange of light 

{M^f^^^) and heavy neutrino (M^"^"^^) must be computed in a nuclear 
model. Eq. (^, however, applies to any intermediate particle. 

At this point we should stress that the main suppression in the mass 
terms comes from the smallness of neutrino masses. In the case of heavy 
neutrino not only from the large values of neutrino masses but the small 
couplings, (7^^^^ for the left handed neutrinos and k and e for the right- 
handed ones. 



2.2. The leptonic left-right interference mechanism (X and 77 terms). 

As we have already mentioned in the presence of right handed currents one 
can have interference between the leptonic currents of opposite chirality. 
This leads to different kinematical functions and two new lepton violating 
parameters A and 77 defined by 

3 

V = evRL , X = KVRL , VRL ^ (U^l'^U^l'^) (8) 

1 

The parameters A and 77 are small not only due to the smallness of the 
parameters k and e but in addition because of the smallness of C^(2i)- 

All the above contributions vanish in the limit in which the neutrino is 
a Dirac particle. 

Many nuclear matrix elements appear in this case, but they are fairly 
well known and they are not going to be reviewed here (see e.g. ^ and 

^, ^, ^ and in our notation ||l^). We only mention that in the case of 
the 77 we have additional contributions coming from the nucleon recoil term 
and the kinematically favored spin antisymmetric term. These dominate 
and lead to values of 77 much smaller than A (l^ . 
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3. The R-parity violating contribution 

In SUSY theories R-parity is defined as 

with B — baryon, L = lepton numbers and s the spin. It is +1 for ordinary 
particles and -1 for their superpartners. R-parity violation has recently 
been seriously considered in SUSY models. 

R-parity violating terms may induce a Majorana neutrino mass and 
may, therefore, lead to Oi^/3/3 decay. But the relevant masses are small, 
at least 10 times smaller |^ than those deduced from the present data. 
The bilinear terms in the superpotential also lead to mixings between the 
neutrinos and neutralinos as well as between the leptons and the charginos, 
leading to lepton violating processes |Q. We are not going, however, to 
consider such effects in this review. 

Here we will be concerned with trilinear couplings in the superpotential 
given by: 

W = \,,km';E',eab + K^k^U^Dleab + K^kU^U^Dl (10) 

where a summation over the flavor indices i,j,k and the isospin indices a,b 
is understood ( Ayfe is antisymmetric in the indices i and j) The last term 
has no bearing in our discussion, but we will assume that it vanishes due 
to some discreet symmetry to avoid too fast proton decay. The A's are 
dimensionless couplings not predicted by the theory. 

In the above notation L,Q are isodoublet and i?'^, D'^ isosinglet chiral 
superfields, i.e they represent both the fermion and the scalar components. 
It has been recognized quite sometime ago that the second term in the 
superpotential could lead to neutrinoless double beta decay Q and 
re-examined quite recently ||l7| . Typical diagrams at the quark level are 
shown in Fig.l. Note that as intermediate states, in addition to the s- 
leptons and s-quarks, one must consider the neutralinos, 4 states which are 
linear combinations of the gauginos and higgsinos, and the colored gluinos 
(supersymmetric partners of the gluons) . Whenever the process is mediated 
by gluons a Fierz transformation is needed to lead to a colorless combina- 
tion. The same thing is necessary whenever the fermion line connects a 
quark to a lepton. As a result one gets at the quark level not only scalar 
(S) and pseudoscalar (P) coupHngs, but tensor (T) couplings as well. This 
must be contrasted to the V and A structure of the traditional mechanisms. 
One, therefore, must consider how to transform these operators from the 
quark to the nucleon level. 

The effective lepton violating parameter, assuming that pion exchange 
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Figure 1. The R-parity violating contribution to Qv^P decay mediated by 
s-fermions and neutralinos (gluinos). 



mode dominates, as the authors of Ref g] claim, is given by 

.3, 



VsusY = (^iii)^^(XPS VPS + m) 



(11) 



with rjps{riT) associated with the scalar and pseudoscalar (tensor) quark 
couplings given by 

VPS = Vx,e + Vx.q + ?7x,/ + Vg + ^V'g , VT = Vx.q - ^x,/ + ^~9-^'g (12) 

They find xps — (2/3), but as we shall see it depends on ratios of nuclear 
matrix elements. For the diagram of Fig. la one finds 



Vx,e 



2 \2 ('*^) ^'^^ 
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(13) 



For the diagram of Fig. lb one finds 



For the diagram of Fig.lc one finds 



% f " ^777^ 2-T2 '^'^^d ed + ^^>^u{ )eu + K^Kfi 16) 

2(Gi^TO^)^ m^^ to^ 



TT 1 TO 

where 



iz! (Gri?TO^yj^ TOg 



mx TOj^ 
4 

(4)//'= E ^..Z^../';;^ (19) 

where e~ f and e~. ?/ are the couphngs of the i^^ neutrahno to the relevant 
fermion-sfermion, which are calculable ( see e.g Ref. |57|). Thus ignor- 
ing the small Yukawa couplings coming via the Higgsinos and taking into 
account only the gauge couplings we find 

^ Z,, + tan6wZ,, 
sinOw 

_ + {tanOw l'i)Zii _ Zu 

where Z\i^Z2i are the coefficients in the expansion of the B,W3 in terms 
of the neutralino mass eigenstates. Note that in this convention some of 
the masses m^^. may be negative. 



4. The effective nucleon current 

As we have mentioned the effective nucleon current in addition to the usual 
V and A terms (P,S,T in SUSY contributions) contains momentum depen- 
dent terms |^ . 

Within the impulse approximation the nuclear current in Eq. (|l|) 
expressed with nucleon fields 5" takes the form 



J! 



2m, 



*, (22) 
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where M is the nucleoli mass, = {p — p')fi is the momentum trans- 
ferred from hadrons to leptons (p and p' are four momenta of neutron 
and proton, respectively) and af"" = (i/2)[7'', 7"]. gv{q^), gM{q^), gA{q^) 
and gp{q'^) are real functions of a Lorenz scalar . The values of these 
form factors in the zero-momentum transfer limit are known as the vector, 
weak-magnetism, axial vector and induced pseudoscalar coupling constants, 
respectively. 

For nuclear structure calculations it is necessary to reduce the nucleon 
current to the non-relativistic form. We shall neglect small energy transfers 
between nucleons in the non-relativistic expansion. Then the form of the 
nucleon current coincides with those in the Breit frame and we arrive at 



A 

J'^{x) = Y.<^9''''j\q^)+9''''Jl{q^mx~?n). fc = l,2,3, (23) 

n=l 

with 

J\q') = 9v{q% = gM{q')^"-^ + gA{q')[S-^^] (24) 

2M q'^ + 

fn is the coordinate of the nth nucleon. 

For the form factors we shall use the following parameterization ||2^ : 

9v{q^) =5^/(1 + 9 VAv)^ gniq^) = i^tp - f^n)gviq^), 

gA{q^)^gA/{l + qV^\? 

where gv = 1, 5a = 1-254, {^x.p - /i„) = 3.70, = 0.71 {GeVf |§ 
and = 1.09 GeV [Q. In previous calculations only one general cut-off 
Ay — Aa ~ 0.85 GeV was used. In this work we take the empirical value 
of A^ deduced from the antineutrino quasielastic reaction I7^p ji^n. A 
larger value of the cut-off A^i is expected to increase slightly the values of 
corresponding nuclear matrix elements. It worth noting that with these 
modifications of the nuclear current one gets a new contribution in the 
neutrino mass mechanism, namely the tensor contribution. The two body 
effective transition operator takes in momentum space the form 

n = T+r+{-gviq'^) + her cr^ - hTSi2) (25) 

where the three terms correspond to Fermi (F) , Gamow- Teller (GT) and 
Tensor (T). One finds that 

S'12 = 3((Ti • gCT2 • 9) - 0-12, (T12 = CTi • (72. (26) 

Note that the tensor operator is defined in momentum space ( q rather 
than f) and there is a change of sign in going to the coordinate space. 

hGT{q') ^ e . 1 , e .2 1 , 2 gUP)e 

g\{q^) ^ i q-^+ml i ^q^ + ml' ^ 3 4m2 ' 
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Table 1. The Fermi, Gamow- Teller and Tensor nuclear matrix elements for 
the light Majorana neutrino exchange of the Oi'/3/3-decay of ^^Ge and ^^°Te 
with ( rows 2 and 4) and without (rows 1 and3) short-range correlations. 



transition 


Gamow- Teller 


Tensor 










AA 


AP 


PP 


AP 


PP 








5.132 


-1.392 


0.302 


-0.243 


0.054 


-2.059 


4.042 


-0.188 




2.797 


-0.790 


0.176 


-0.246 


0.055 


-1.261 


2.183 


-0.190 




4.158 


-1.173 


0.258 


-0.329 


0.074 


-1.837 


3.243 


-0.255 




1.841 


-0.578 


0.134 


-0.333 


0.075 


-1.033 


1.397 


-0.258 



Table 2. Nuclear matrix elements for the light and heavy Majorana neu- 
trino exchange modes of the 0:^/3/3-decay for the nuclei studied in this work 
calculated within the renormalized pn-QRPA. 

{30)t)iy ~ decay : 0+ — > 0+ transition 
M. E. s^si rwrpo rre^rf ns^i T^Wi rasxi istr^ 









liglit Maj 


orana 


neutrino (1= 


:ligllt) 








Kv 

Ka 
Kp 
Kp 


0.80 


0.74 


0.45 


0.82 


0.50 


0.75 


0.66 


0.32 


1.14 


2.80 


2.66 


1.54 


3.30 


2.08 


2.21 


1.84 


0.70 


3.37 


0.23 


0.22 


0.15 


0.26 


0.15 


0.24 


0.21 


0.11 


0.35 


-1.04 


-0.98 


-0.65 


-1.17 


-0.69 


-1.04 


-0.91 


-0.48 


-1.53 


Km.> 


2.80 


2.64 


1.49 


3.21 


2.05 


2.17 


1.80 


0.66 


3.33 








heavy Majorana 


neutrino (1= 


heavy) 








Kv 

Ka 
Kp 
Kp 


23.9 


22.0 


16.1 


28.3 


17.2 


25.8 


23.4 


13.9 


39.4 


-55.4 


-51.6 


-38.1 


-67.3 


-39.8 


-60.4 


-54.5 


-31.3 


-92.0 


106. 


98.3 


68.4 


123. 


74.0 


111. 


100. 


58.3 


167. 


13.0 


12.0 


9.3 


16.1 


9.1 


14.9 


13.6 


7.9 


23.0 


-55.1 


-50.7 


-41.1 


-70.1 


-39.0 


-64.9 


-59.4 


-34.8 


-101. 




32.6 


30.0 


14.7 


29.7 


21.5 


26.6 


23.1 


14.1 


35.6 



hriq^) ^ . 2 _ 1 ^ 1 gW^W^ .g?^ 

gliq^) ^3q^+ml 3^q^ + ml'^ ^3 Am^ ' ^ ' 

The exact results will depend on the details of the nuclear model, since 
the new operators have different momentum (radial) dependence than the 
traditional ones and the tensor component is entirely new. We can get a 
crude idea of what is happening by taking the above average momentum 
(g)=100 McV/c. Then we find that the GT ME is reduced by 22%. Then 
assuming that T matrix element is about half the GT one, we find that 
the total reduction is 28%. This is in perfect agreement with the exact 
results for the A~76 system, 29%, but a bit smaller than the 38% obtained 
for the A=130 system. We will now summarize the results obtained with 
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the above modifications of the nucleon current for hght neutrino (see Eq. 

f) for the two representative Oi^/3/3-decay nuclei "^^Ge and ^^^Te in Table 
The details of our calculations will be given elsewhere One notices 
significant additional contributions to GT (AP and PP) and tensor (AA 
and PP) nuclear matrix elements coming from higher order nucleon current 
terms. AP and PP originate from the second (first) and third (second) 
terms in /igt(^t) of Eq. (p7|). 

By glancing at the Table ^ we also see that, with proper treatment 
of short-range two- nucleon correlations (see e.g Vergados |^), all matrix 
elements are strongly suppressed. The effect is even stronger in the case 
of heavy intermediate particles. Detailed results [ p3| for various nuclei are 
presented in Table ||. 



5. Extraction of the lepton violating parameters 

The limits deduced for the Icpton-number violating parameters depend on 
the values of nuclear matrix element, of the kinematical factor and of the 
current experimental limit for a given isotope [see Eq. (^]. 

5.1. Traditional lepton violating parameters 

Even, though, we expect the nuclear matrix elements entering the light 
neutrino mass mechanism to be decreased by about 30% , independently 
of the nuclear model, we will stick to the calculations as reported. Thus 
the present best experimental limits ^ |^ can be converted to upper 
limits on < > and r/j^ . 

The thus obtained results are given in Table The references of Ta- 
ble H are defined as follows: Exl=Heidelberg-Moscow Collaboration |61], 
Ex2=Elliott et al [6^, E x3^Kawashima et al jel, Ex4=Ejiri et al |4], 
Ex5=Davenich et al [|65[, Ex6=Bernatovicz et al |6ra, Ex7=Alessandrello 
et al^, Ex8=De Silva et al Ex9=Busto et al ^ Thus, the most 
restrictive limits are as follows: 



< nil, > 



best 



< 0.62 eV 



best 



< 1.0 X 10" 



(28) 



(see Ref.|23|, ||6l|) By assuming < m^, >=< and = 77^^*=^*. 

(1) we calculated half-lives of the Oi//3/3-decay T^J^'^^iK m„ >6e<*t)^ 
rjnexp-Oi' ^^_^best^ for nuclcar systems of interest using specific mechanisms 
with the "best" parameters. The thus obtained results are given in Ta- 
ble m The references of Table ^ are defined as follows: Exl=Heidelberg- 
Moscow Collaboration] 6l| , Ex2=Elliott et al ^2\, Ex3=Kawashima et al 
[p3 |, Ex4=Ejiri et al 1^4 1, Ex5=Davenich et al |65|, Ex6=Bernatovicz et al 
|§, Ex7=Alessandrello et al Ex8=Busto et al Ex9=De Silva et 
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Table 3. The present state of the Majorana neutrino mass searches in 
/3/3-decay experiments. T^'^J^~°'' (present) is the best presently available 
lower limit on the half-life of the Oz//3/3-decay for a given isotope. The 
corresponding upper limits on lepton number non-conserving parameters 
< nil, > and 77 jv are presented. For the definition of the references and 



best" see text. 



Nucleus 


76 Ge 


«2Se 


96 Zr 


100 Afo 


"6cd 


T-^-"" (present) [y] 
Ref. 

<mr, > [oV] 

T^jr [y] 


1.1 X 10^5 

[Exl] 
0.62 

1.1 X 10^5 

1.0 X 10~7 

1.1 X 10^5 


2.7 X 10^2 
[Ex2] 

6.3 

2.8 X 10^" 
1.1 X 10^6 

2.9 X 10^4 


3.9 X 10^9 

[Ex3] 

203. 

4.2 X lO^** 
4.0 X 10-'' 
5.8 X lO^-* 


5.2 X 10^2 

[Ex4] 

2.9 

1.2 X lO^-* 
6.2 X 10-7 
1.8 X lO^-* 


2.9 X 10^2 

[Ex5] 

5.9 

2.6 X 10^4 

1.1 X 10-6 

3.2 X 10^* 


Nucleus 






i«6Xe 






T-^''-"" (present) [y] 
Ref. 

<m^ > [eV] 


7.7 X 10^* 

[Ex6] 

1.8 


8.2 X 10^1 

[Ex7] 

13. 


4.2 X 10^3 

[Ex8] 

4.9 


1.2 X 10^1 

[Ex9] 

8.5 





rpexp—Qu 

1/2 

(< >6esi) [y] 6.6 X 10^5 3.8 X 10^4 2.7 x lO^"' 2.3 x 10^3 

r]N 2.9 X 10-7 2.0 X 10-6 4.5 x 10-7 1.6 x 10-6 

T^jr^'ly^i^N"') M 5.9x1025 3.1x1024 7.9x1024 2.7x1023 



al |68|, Since the quantities < rui, >, ?7„ depend only on particle theory 
parameters these quantities indicate the experimental half-life limit for a 
given isotope, which the relevant experiments should reach in order to ex- 
tract the best present bound on the corresponding lepton number violating 
parameter from their data. Some of them have a long way to go to reach 
the Ge target limit. 

A summary involving most of the available nuclear matrix elements and 
taking into account what, at present, is a good guess as canonical values 
of the lepton violating parameters is provided in Tables ^ and The 
references in these tables are defined as follows: R— Retamosa et al p^ , 
H^Haxton et al [|o|, El=Engel et al E2=Engel [||, S=Suhonen 
et al^, M=Muto @, T=Tomoda et 0, Pl^Pantis et a? 
P2=Pantis et al (p-n pairing), Sl=Simkovic et al (and private com- 
munication), F=Faessler et al Q, Q P==Present calculation (see Simkovic 
et al ||2^ for the nuclear Matrix elements). Notice in particular that the 
present calculation, marked P in the table, involves not only renormal- 
ized QRPA but takes into account the corrections in the hadronic 

current discussed above (see table 
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Table 4. The lifetimes predicted for 0+ — > 0+ 0^(3(3- decay in various mech- 
anisms (light neutrino, heavy neutrino, A and ry terms and SUSY contribu- 
tion) for suitable input of lepton violating parameters and available nuclear 
calculations. For the definitions of the references see text. 

— decay : 0+ 0+ transition 









^1/2 


"'^'"^i{m,),(X),{v),{vN),{vsuSY}) 


















iibc^ i^»Te ""Te 






Rcf. 


10^4 


W'^ 


10^" 


10 10 


10 10 


10 


10^-* 








(mv) = 


= leV, (A) 


= 0, = 0, (r?jv> 


= 0, (risusv) = 








R 


12.8 


34.8 


4.80 








24.2 




H 


6.34 


3.36 


1.16 




0.80 


0.32 






El 




4.60 


1.84 




0.90 


0.48 






E2 




28.0 


11.2 




3.00 


1.32 


6.60 




S 




8.12 


2.86 




3.60 


1.66 






M 




4.66 


1.20 


2.54 


1.54 


0.98 


4.42 


6.74 


T 




4.32 


1.22 


0.52 


1.96 


1.08 


2.80 


8.90 


PI 


5.00 


7.20 


3.00 


1.22 7.80 


9.40 3.80 


1.72 


6.60 




P2 


56.0 


36.0 


5.60 


54.0 


9.80 30.0 


4.20 


5.60 




SI 




17.9 




0.50 


1.44 2.18 




17.5 




P 




4.22 

{mv) = 


1.08 

0, (A) = 


1.61 0.46 
10-«,(r?> =0, (r?jv; 


0.99 2.53 

> = 0, (vsusy) = 


1.46 


10.1 


8.78 


R 


7.45 


50.2 


3.25 








22.2 




S 




7.75 


1.14 




14.8 


0.89 






M 




7.35 


0.99 


0.95 


13.5 


0.95 


4.90 


3.73 


T 




8.02 


1.07 


0.55 


21.1 


1.18 


3.47 


6.71 


PI 


2.71 


8.90 


2.08 


0.94 30.6 


39.1 22.7 


1.34 


2.73 




P2 


27.9 


41.2 


4.39 


27.7 10.3 


10.8 165 


2.22 


4.42 





5.2. R-parity induced lepton violating parameters 

In this section we will elaborate a bit further on the R-parity violating 
parameters. We will consider the pionic contribution (^. We will first 
attempt to evaluate the relevant amplitude using harmonic oscillator wave 
functions, but adjusting the parameters to fit related experiments. 

Let us begin with the second process of Eq. (|l]) . This process involves a 
direct term and an exchange term. The direct term is nothing but a decay 
of the pion into two leptons with a simultaneous change of a neutron to a 
proton by the relevant nucleon current, which in this case can only be of 
the of the PS type. The tensor contribution cannot lead to a pseudoscalar 
coupling at the nucleon level, which is needed to be coupled to the usual 
pion nucleon coupling in the other vertex to get the relevant operator for 
a 0"*" — > 0"^ decay. Thus the amplitude involving the meson is related to 
TT, ^ decay. 

Ai.^{direct) = 4 ai-^ J^^'^ ml exp{-^^-^^ (29) 
(2 toat) 6 
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Table 5. The previous table continued. 





{tTIi/ ) 


— n l\\ 


— w/ — J^*-* ) WiV 


) — 0) {vsusy) 


— n 

— u 






R 


6.42 


27.2 


6.24 






22.2 








OU. 1 


11 1 


10 7 








IVl 




i .oo 


n QQ n Qf^ 


iO.O 


u.yo 


A QO 


O. i o 


T 




2.25 


0.65 0.28 


0.67 


0.44 


1.21 


3.39 


PI 


15.11 


3.10 


6.51 1.48 3.44 


19.2 1.20 


0.62 


1.23 




P2 


43.2 


22.8 


5.16 7.95 102 


83.2 1.90 


1.05 


0.96 








= 0, (A) 


= 0, (»;> = 0, im) = 


10"'^, {risusv) 


= 






PI 


4.95 


0.25 


3.35 67.1 4.70 


23.5 0.78 


3.03 


1.42 




P2 


124 


0.59 


7.23 671 1.47 


33.6 1.27 


1.31 


1.01 




P 




15.4 


4.10 8.10 0.97 


8.40 8.51 


4.54 


3.94 


40.6 




{rriv) = 


= 0, (A) = 


= 0, {ri) = 0, ( etajv) = 


= 0, (tisusy) = 


io-« 






F 




3.3 


0.86 0.71 0.30 


0.85 0.93 


0.45 


1.2 


3.2 


P 




4.5 


1.0 1.4 0.59 


1.4 1.5 


0.71 


2.0 


5.2 



with 

ai, = (27r)3 ^ x(0) (30) 

The exchange contribution, in which the produced up quark of the 
meson is not produced from the ^\^acuum.^^ but it comes from the initial 
nucleon, is a bit more complicated. The harmonic oscillator quark model, 
however can be used to get its relative magnitude (including the sign) with 
respect to the direct term. This way we find 

Ai^{exchange) = -3 ai^[-—-^——ml exp{ — ] (31) 

Thus the effective two-body transition operator in momentum space at the 
nucleon level becomes: 



PS _ r ^i-g^2.g , (qb)^ 

^^'^ ~ ''^^ \2mNf 6 \-' + ml ^^^^ 



with CiTT = 5r Oiiv 

ci. = (2^)' ^ 9r X(0) (33) 

where Qr = 13.5 is the pion nucleon coupling and niq is the constituent 
quark mass. Note the presence of the exponential form factor in the har- 
monic oscillator model, which has been ignored in other treatments. We 
see that, in going from the quark to the nucleon level, the factor of three 
coming from the mass gain is lost due to the momentum being reduced by 
a factor of three. The quantity x(0) is essentially the meson wave function 
at the origin given by: 

X(0) = rn-J'' m (34) 
16 



The quantity x(0) can be obtained from the n fi,!/ decay via the expres- 



i^i(%^m^m^(l-4)'x^(0) (35) 



sion 

1 ^,Gf.2 2 
— = — I I m Ti, , , 

ml 

From the measured hfetime r = 2.6 x 10^* we obtain x(0) = 0.46 

The first process of Eq. |l| is easier to handle. Now both the PS and T 
terms contribute. We thus get 

3 1 

^27^(7") = - (527r rn^ ,A2APS) = - a2-K , (36) 

52. = 4 (2^)3 x\0) (37) 

T _ ai.ga2.g ml pg _ 2 ^ 

(2 mi^Y {q-^ + m^j 3 

- 4 (2^)3 x'(0) , C2. = 4 (27r)3 .g^ ^^(O) (39) 

Using the above value of x(0) and — 13.5 we get cx-n — 109 and 
C2. = 198 which are in good agreement with the values 132.4 and 170.3 
respectively obtained by Faessler et al JlSf . 

It is now customary, but it can be avoided to go to coordinate 
space and express the nuclear matrix elements in the same scale with the 
standard matrix elements involving only nuclcons. Thus we get 

ME^ = (^)2 au. ^ \M^^ + M|-] (40) 

Where the two above matrix elements are the usual GT and T matrix 
elements with the additional radial dependence given by 

i^^^ =e-^ , F^'^ = (3 + 32; + 2;2) ^ (41) 

Fl\ = (x - 2)e-^ , Fl"" (1 + x) e-^ (42) 



aiTT = - ci^p,a2. = C27r p , P=7X72 ( ) ( ) (43) 

48/^ TOp mA 

In the above formulas we have tried to stick to the definition of r\suSY 
given above (see |l^), but since the tensor (at the quark level) does not 
contribute to the Itt diagram the "effective" nuclear matrix element is not 
the sum of the two matrix elements of Eq. (|4^), but only ME2, and xps 
depends on the nuclear matrix elements, i.e. 

2 AIE 

ME^ff = ME2 , = 3 (4^ + 1) (44) 

17 



There is no difference, of course, between the two expressions if ME2 is 
dominant, as is actuaUy the case). 

Before proceeding further we should remark that for the experimen- 
taUy derived harmonic oscillator parameter for the tt meson, b — 1.8 /, 
a27r s dominant and xps approaches the value of 2/3. In fact we find 
aiTT = — 1.2 X 10^^ and ai^r = 0.15 which are in good agreement 
with the values —4.4 x 10^^ and 0.20 respectively obtained by Faessler et 
al [|l^. Furthermore from the nuclear matrix elements of Ref. one 
can see that the M"^^ is favored, since, among other things, its tensor and 
Gamow- Teller components are the same magnitude and sign (in the Itt 
mode they are opposite). Thus nuclear physics also favors the 2tt mode. 

With the above ingredients and using the nuclear matrix elements of 
we can extract from the data values of rjsusY- 

Then one can use these values of rjsuSY in order to extract values for 
the R-parity violating parameters X'lu. 

As we have already mentioned one must start with 5 parameters in the 
allowed SUSY parameter space and solve the RGE equations to obtain the 
values of the needed parameters at low energies ^ Q . For our purposes 
is adequate to utilize typical parameters, which have already appeared in 
the literature |7^. One then finds 

X'lii — Cj^o{risusYy^^ [neutralinos only) (45) 

X'lii^ CgivsusvY^'^ [gluino only) (46) 
When both neutralinos and gluinos are included we write 

Kii ^C^o~^{risusY)"^ (47) 

The values of these coefficients are given in Tab. 6 for the nine SUSY 
models mentioned above. From Table 6 we see that there is quite spread 
in the quantities Cj^o,C'g and Cy^O g, depending on the SUSY parameter 
space. We will see that this is the largest uncertainty in estimating the 
SUSY contribution to Oiz/Sf] decay. In all of these cases the intermediate 
selectron-neutralino mechanism appears to be the most dominant. The most 
favorable situation occurs in the case # 7 of Table 6. And this what we 
will consider in extracting the limits on X'm Combining the above values of 
the couplings ctfe^, k=l,2, with the corresponding nuclear matrix elements 
of of Faessler et al |l|l) (F) and the two nucleon ME of Wodecki et al @ 
(W) we obtain the limits listed as Pr in Table 7. 

Thus the most stringent limit is obtained from the ''^Ge data and is 

A'ni < 6.0 X 10"^ {for case #7 ) (48) 

The above quantities are assumed positive. If not, the absolute value is 
understood. 
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Table 6. A sample of relevant parameters obtained by some choices in 
the allowed SUSY parameter space. It clear that in all cases the neutralino 
mediated mechanism is dominant (for definitions see text) . The parameters 
C shown have been multiplied by 10~^ 



input Kane et al (# 1-3) Ramond et al {# 4-9) 





# 1 


#2 


#3 


#4 


#5 


#6 


#7 


#8 


#9 


tan/B 


10. 


1.5 


5.0 


5.4 


2.7 


2.7 


5.2 


2.6 


6.3 


mo 


124. 


26 


96 


83 


124 


58 


34 


34 


50 


mo 


237. 


65 


173 


150 


204 


108 


66 


74 


92 


mo 

^3 


455. 


219 


310 


391 


445 


336 


170 


191 


208 


mo 


471. 


263 


342 


409 


472 


361 


208 


236 


244 




328. 


124 


211 


426 


472 


310 


90 


94 


109 


ma^ 


700. 


283 


570 


590 


664 


449 


251 


275 


319 


my 
m-g 


676. 


276 


550 


577 


638 


441 


246 


268 


310 


718. 


292 


610 


483 


706 


371 


280 


304 


350 


C^a X 10"^ 


3.3 


0.023 


0.46 


5.9 


14 


1.4 


0.0068 


0.0089 


0.019 


C-g X 10-3 


14 


1.6 


54 


56 


110 


13 


0.97 


1.5 


3.1 


C.o g X 10-3 


3.2 


0.023 


0.45 


5.3 


12 


1.3 


0.0068 


0.0089 


0.019 



Table 7. The limits for rjsusY and X'm obtained : a) For the pion mecha- 
nism using the values of aiT^ and a2Tr computed in this work and the nuclear 
ME of Faessler et al (F) b) Using the nuclear ME of the two nucleon mode 
of Ref. Wodecki et al (W). In extracting the values of A'^^^ we used the 
SUSY data of #7 of Table 6. The experimental lifetimes employed are 
those of Table 4. 



Pion mode Only nucleons 



(A,Z) 


risuSviPr) 




VsusriF) 




VSUSviPr) 






8.4 X 10-^ 


6.0 X 10' 


4 


5.5 X 10-^ 


4.8 X 10- 


'4 


2.6 X 10-* 


1.1 X 10-3 


100 Mo 


3.2 X 10-* 


1.8 X 10- 


'3 


2.4 X 10-* 


1.1 X 10- 


-3 


1.1 X lO-'^ 


2.2 X 10-3 


116 


7.6 X 10-* 


1.8 X 10- 


'3 


5.4 X 10-* 


1.5 X 10- 


^3 


2.6 X 10-'^ 


3.3 X 10-3 


128 


1.6 X 10-* 


8.1 X 10- 


-4 


1.1 X 10-* 


6.8 X 10- 


-4 


5.6 X 10-* 


1.6 X 10-3 


iSOTe 


1.0 X 10-^ 


1.8 X 10- 


-3 


5.5 X 10-* 


1.5 X 10- 


^3 


1.3 X lO-'^ 


6.5 X 10-3 


l36Xe 


2.4 X 10-* 


9.4 X 10- 


^4 


1.7 X 10-* 


7.8 X 10- 


^4 


8.7 X lO-'^ 


2.2 X 10-3 


isoATrf 


7.3 X 10-* 


2.3 X 10- 


-3 


5.2 X 10-* 


1.4 X 10- 


-3 


2.4 X lO-'^ 


3.0 X 10-3 



6. Conclusions 

We have seen that Oi>f3P decay pops up in almost any fashionable particle 
model. Thus it can set useful limits not only on the light neutrino mass 
(28), but in addition on other lepton violating parameters like < 77 at > of 
(28) or the parameters A and 77 (see sect. 2.2). Finally we mention again 
the limit extracted on the R-parity violating parameter (^). A set of 
limits, for our choice of nuclear matrix elements, derived from the various 
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Table 8. Summary of the results presented in this work. 





{irijy) 


(A) 


in) 


(vn) 


{'nsusv) 


A'ni 


(A,Z) 


eV 


lOx"" 


10x-» 


10x-» 


10x-» 


lOx-* 




Pr 


PI 


PI 


Pr 


Pr 


Pr 




0.27 


0.56 


0.32 


0.44 


0.31 


4.0 




2.9 


26 


8.8 


6.2 


3.2 


18 


i«cd 


5.9 


37 


26 


11 


7.6 


18 


i28Te 


1.8 


5.6 


1.3 


2.9 


1.6 


8.1 


i30Te 


13 


7.6 


5.2 


20 


10 


18 




49 


2.1 


1.4 


4 5 


2.4 


9.4 


150 Afd 


8.5 


5.6 


5.3 


16 


7.3 
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nuclear targets is given in Table 8. For ^^Ge our results are different from 
the ones given above, since we have used here the unpublished new limit 
of the Heidelberg- Moscow experiment T1/2 > 5.7 x 10^^ yr. 

We see that limits are quite stringent, but they, of course, have uncer- 
tainties in them. They come from nuclear physics, especially for the short 
ranged operators or from particle physics, as, e.g., in the case of super- 
symmetry. It is however evident that in the extraction of X'ln the main 
uncertainty comes from the parameters of supersymmetry. After all X'm 
depends on the inverse fourth root of the lifetime and the inverse square 
root of the nuclear matrix elements. On the other hand it depends in the 
second power of the masses of the mediating SUSY scalars. This is not 
true of the models considered here but it is found in other calculations as 
well 0. 

It is clear that during the last year the interest of most people is be- 
ing focused on the light neutrino mass mechanism. This due to the ex- 
perimental indications for neutrino oscillations of solar (Homestake | |75[ ], 
Kamiokand c [|76| , Gallex |Q and SAGE ||7^), atmospheric (Kamiokande 
dll, 1MB fsj and Soudan Super-Kamiokande |||) and terrestrial 
(LSND experiment [Q) experiments. 

One can use the constraints imposed by the results of neutrino oscil- 
lation experiments on < Ato^ >. These experiments, of course, cannot 
predict the scale of the masses or the Majorana phases. The predictions 
differ from each other due the different input and structure of the neutrino 
mixing matrix and assumptions. Bilenky et al and others [ |86[ | have 
shown that under quite reasonable assumptions in a general scheme with 
three light Majorana neutrinos and mass hierarchy | < rrii, > \ is smaller 
than 10~^ eV. In another study outlined in Ref.|^5| the authors end up 
with I < nil, > I « 0.14 eV. Thus one can see that, the current limit on 
< rrii, > in (28) is quite a bit higher than the neutrino oscillation data. 

There is a new experimental proposal for measurement of the Qv(3(3- 
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decay of ^^Ge, which intents to use 1 ton (in an extended version 10 tons) 
of enriched ^^Ge and to reach the the half-hfe hmit T^l'^'^'^ > 5.8 x 10^^ 

and T'^J~'^^^ > 6.4 x 10^^ after one and 10 years of measurements, respec- 
tively. From these half-hfe values one can deduce [see Eq. d) and Table 
^ the possible future limits on the effective light neutrino mass 2.7 x 10~^ 
eV and 8.1 x 10"'^ eV, respectively. From the comparison with limits advo- 
cated by the neutrino oscillation phenomenology we conclude that GENIUS 
experiment |6^, would be able to measure this lepton number violating 
process, provided, of course, that the neutrinos are Majorana particles. 

We must emphasize that the plethora of other 0i/f3(3-deca.y mechanisms 
predicted by GUT's and SUSY do not diminish the importance of this 
reaction in settling the outstanding neutrino properties. One can show 
that the presence of these exotic mechanisms implies that the neutrinos 
are massive Majorana particles, even if the mass mechanism is not the 
dominant one |8). 

Thus one can say with certainty that the experimental detection of the 
Oj^/3/?-decay process would be a major achievement with important impli- 
cations on the field of particle and nuclear physics as well as on cosmology. 
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